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ABSTRACT 

The present paper deals with the speed estimation of the induction motor using observer with 
Model Reference Adaptive System (MRAS). The first part of the paper includes the mathematical 
description of the observer for the speed estimation of the induction motor. The second part 
describes the simulation results which are shown for different changes of the induction motor speed 
which confirm high dynamic properties of the induction motor drive with sensorless control. The 
third part describes experimental implementation using DSP controller. An adaptive state observer, 
MRAS is tested on two different schemes to observe rotor speed. The stability analysis and 
continuous and discrete models of each scheme are investigated. The high performance of these 
schemes is shown in simulations and experimental results. 

Keywords: Model Reference Adaptive System (MRAS), Sensorless Control, reactive power 

I. INTRODUCTION 

Sensorless control techniques for IM drives have been widely investigated over the last two 
decades. The great advantages offered by sensorless control including compactness and robustness 
make it attractive for many industrial applications specially those operating in hostile environments. 
Among several strategies proposed for sensorless IM drives, model reference adaptive systems 
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(MRAS) are the most popular schemes employed due to their simple implementation and smaller 
computational effort. However, these schemes usually fail to provide a satisfactory response at low 
stator frequency. Much research interest has been devoted to improve the performance of MRAS- 
based sensorless schemes in this region of operation. Adaptive system can be defined as a system 
that consists of a primary feedback that takes care of process signal variations and a secondary 
feedback that deals with process state changes. In this definition, the primary feedback is used as in 
non-adaptive control, and the secondary feedback makes the system adaptive. From this definition it 
is clear that process state variations give rise to adaptation of the system. The aim of reacting to state 
changes is to attempt to maintain a high system performance, even if the process states are unknown 
or varying [2] . In this paper model reference adaptive system is applied to induction motor drive as a 
state observer. 

II. MRAS BASED SPEED ESTIMATOR SCHEME USING SPACE VECTOR 

The model reference adaptive system is one of the most successful adaptive control 
techniques applied to motor control and parameter estimation [2-3], [4-6]. In a MRAS system, some 
state variables, x d ,x q (e.g. back e.m.f components (e md , e mq ) reactive power components (q md ,q mq ) , 
rotor flux components (*T rd , 4\ q ) etc.) of the induction machine, which can be obtained as sensed 
variables such as stator voltage and currents, are estimated in reference model. They are then 
compared with state-variables x~ d and x' q estimated by using adaptive model. The difference between 
these state-variables is then used in adaptation mechanism, which outputs the estimated value of the 
rotor speed (w r ) and adjusts the adaptive model until satisfactory performance is achieved. Fig. 6. 2 
shows MRAS based speed estimator scheme using space vector. Two of the schemes will be 
discussed in the following sections: reactive power and back e.m.f errors are used as speed tuning 
signals. In these expressions i s _ and i u denote the stator voltage and stator current space vectors 
respectively in the stationary frame e m denotes the back emf space vector also in stationary reference 
frame as Aem = e m - e~ m . The symbol A denotes the quantities estimated by the adaptive model. 
Recent MRAS algorithms mentioned in this work avoid both pure integrators and low-pass filters. 
Reactive power scheme described below is robust to both stator and rotor resistance variations, and 
can even be applied at very low speeds (but not zero speed). Both of the observers (reactive power 
and back emf schemes) described below use monitored stator currents and stator voltages. In a 
voltage-source inverted-fed drive, however, it is not necessary to monitor the dc link voltage and the 
stator voltages since the latter can be reconstructed by using the inverter switching states. 




>- 

► 

► 



Reference 

Model 




-> 

*■ 

► 



v = [e d ,e q ] r 



Adaptive 

Model 



X 



q 



w 



Adaptation 

Mechanism 



Fig.l MRAS based speed estimator scheme 
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Fig.2 Coordinates in stationary reference frame 
Equations for an induction motor in the stationary frame can be expressed as: 
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where w r is a vector whose magnitude is rotor electrical angular velocity, and whose direction is 
determined according to right hand system of coordinates as shown in Fig. 2. “®” denotes the cross 
product of vectors respectively. From Eqn.l, e m and structure of MRAS can be derived as follows: 
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We can write the above equations for the direct and quadrature-axis back emf in the following form: 
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III. REACTIVE POWER MRAS SCHEME 



(3) 

(4) 



First let us define a new quantity q m as the cross product of the counter EMF vector e m and 
the stator current 

vector i s . That is, c| = ' - ^ ^ ... ( 5 ) 



q m is a vector, whose direction is shown in Fig.6.3, and whose magnitude q m represents the 
instantaneous reactive power maintaining the magnetizing current. Substituting the Eqns (3-4) for e m 
in Eqn (5) noting that i s ® i s = 0, we have 
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Using Eqn (6) as the reference model and the adjustable model, respectively, an MRAS 
system can be drawn as in Fig. 3, where proportional and integral (PI) operations are utilized as the 
adaptation mechanism. From Eqn (6)), it is evident that the speed estimation system of Fig. 3 is 
completely robust to the stator resistance, besides requiring no integral calculation. 



Reference Model 




Fig.3 System structure of rotor speed observer 



Notice that the representation of complex number is defined for the stator voltages and currents in 
the stationary reference frame i.e. v s = v s d +j v sq and is' = i s d +j i sq 

III.l Reference Model Continuous Time Representation 



The back emf of the induction motor can be expressed in the stationary frame as follows: 
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The reactive power of the induction motor can be computed from cross product of stator currents and 
back emf vectors as follows: 
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As a result the reactive power shown can further be derived as 
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III.2 Adaptive Model Continuous Time Representation 



The estimated back emf computed in the adaptive model can be expressed as follows: 



L m di t m / A . . • \ 

e mq = L ^ -^-(-TrWr'md-^q+lsqJ 
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where, Tr = Lr/Rr 

T = is the rotor time constant, i m d, i m q are computed from the following equations: 
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Once the estimated back emf computed byEqns (17-18), the estimated reactive power can be 
computed as follows: 



*-1111 *s X e m ^sq^md ( 14 ) 

generated by adaptive model matches that generated by reference model. The speed tuning signal is 
the error of reactive power that can be expressed as follows: 



e A e = i s x(e m — e m ) = q m -q m 

If the MRAS successfully maintains nearly zero error, and if the same value of Tr is used in 
the RAS adjustable models and in the function block for calculating w s u p , then we have the 
following relations :w e = w" e and T r w s u p = T" r w' s | ip where variables without “ A ” are actual values, 
and ones with “ A ” represent the corresponding values used in the MRAS vector control systems. 
Thus, if T r ± T% , then w s i ip ^ w' slip , but w 0 = w" 0 , which is used for orienting the stator current 
vector. Therefore, complete field-orientation can be achieved even if the value of T r is quite wrong. 
The error in the value of T r , however, produces an error in the speed feedback, thus affecting the 
accuracy of the speed control as follows: 

Sw = w A r - W r = { 1- T r /TT } Wsiip (16) 

This also holds for the previous MRAS scheme. However, the accuracy of the speed 
estimation system discussed depends on the transient stator inductance and also referred magnetizing 
inductance. The latter quantity is not too problematic, since it does not change with temperature 

III.3 Discrete time representation for DSP based system implementation 

For implementation on DSP based system, the differential equations need to be transformed 
to difference equations. Due to high sampling frequency compared to bandwidth of the system, the 
numerical integration, such as forward, backward, or trapezoidal rules, can be adopted [6]. 

(i) Reference Model 

According to Eqn (10) reference model reactive power is given as : 
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= i v — oL,^ J 

Using backward approximation: 
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And this equation can be further simplified as: 
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where T is the sampling time. 

(ii) Adaptive Model 

According to (14), reactive power in adaptive model is derived as : 
ci = i x e„ = i .e — i e „ . 

• in s in sd mq sq md 



( 18 ) 



whose discrete-time representation is: 
q m ( k ) = i S d( k )e mq (k)-i sq (k)e md (k) 

In order to compute e~ m d (k) and e'mq (k) consider their continuous time representations 
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which have the discrete-time representations as; 
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and i m d(k), i mq (k) can be solved by using trapezoidal integration method which yields continuous 
time representation 
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(iii) Per unit, discrete time representation 

For the sake of generality, the per unit concept is used in all equations. However, for the simplicity 
the same variables are also used in the per unit representations. 

Reference Model 

Dividing (6.34) by base power of Vb lb, then its per unit representation is as follows: 

q„ Ck) = i*,(k)v !q (k) — >^(k)v rf (k) — 

K , < U _ I )i (k> — < k _ 1>) 



Rearranging (25) to have the one in (26), 

q m ( k > = i sd( k )(v sq (k)-K 1 i sq (k-l))-i sq (k)(v sd (k) + K 1 i sd (k-l))pu 

( 26 ) 



Where , Ki= (a L s Ib)/(T Vb) , Vb is base voltage, and R is base current 

Adaptive Model 

Dividing Eqn (21) by base voltage Vb, then yields 

6 m a(k) = K,t-K 3 w r (k)i mq (k)-i md (k) + i^(k)}pu e mq (k) = K,(-K 3 w t (k)i md (k)-i mil (k) + i s<l (k))pu 

and (27) 

L 2 I. L w. 

where K 2 = ^ ™ v . IC 3 = T r w b = — and w b = 27cf b 

is base electrical angular velocity. Similarly, dividing Eqns (28) and (29) by base current R , then 
yields 
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where, 

IV SIMULATION of the MRAS SCHEME 



(29) 

(30) 



The simulations of reactive power schemes is implemented to confirm the theoretical results 
using Matlab/Simulink. In these simulations the voltage and current outputs of induction machine 
model are used as the inputs of MRAS schemes. Two independent observers are configured to 
estimate the components of back emf and reactive power. The error between the outputs of the two 
observers is then used to derive a suitable adaptation mechanism which generates the estimated 
speed for the adaptive model as shown in Fig. 4. In Fig.4 the adaptive model is configured according 
to equation (17) and the reference model is configured according to equation (18). Fig. 5 and Fig. 6 
shows simulation results for four-quadrant speed tracking performance of the reactive power MRAS 
scheme with 20 hp and 5 hp induction motors. The simulations here are implemented under no-load 
conditions 
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Fig-4 The Simulink model of reactive power MRAS scheme 

In Fig. 5 (b) and Fig.6.(b), more accurate speed estimation and speed error (difference 
between the actual speed and estimated speed) are shown. No-load performance of the speed 
observer is very high as seen in these figures even at very low speed range. Since there does not exist 
any immediate transient speed change due to mechanical loading, speed outputs obtained in the 
simulation are very smooth with negligible speed errors. 
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(a) Estimated speed produced due to J and i sq (b) Estimated speed, speed error 

Fig-5 Four-quadrant speed reversal of 5 hp induction motor at no_load up to rated speed 
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(a)Estimated speed, produced torque due to J and i qs 



20 hp ind. motor four quadrant speed tracking 




(b) Estimated speed, speed error 



Fig.6 Four-quadrant speed reversal of 20 hp induction motor at no_load up to rated speed. 



In Fig. 7 and Fig. 8, the generated electromagnetic torque and q axis stator current are shown 
with estimated speed for both 5 hp and 20 hp induction motors at the same loading conditions. 



5 hp ind. motor speed estimation %100 rated torque %100 rated speed 







Time(sec) 

(a) 100 % rated torque & 100 % rated speed 



5 hp ind. motor speed estimation %10 rated torque %100 rated speed 




(b) %10 rated torque & %100 rated speed 



5 hp ind. motor speed estimation %10 rated torque %10 rated speed 



(c)%10 rated torque & %10 rated speed 

Fig.7. 5hp induction motor estimated speed speed using reactive power MRAS 
scheme, applied torque and stator q-axis current 
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20 hp ind. motor speed estimation %100 rated torque %I00 rated 




20 hp ind. motor speed estimation %10 rated torque %100 rated speed 




(a)%100 rated torque, %100 rated speed (b) %10 rated torque, %100 rated speed 



20 hp ind. motor speed estimation %tO rated torque %tO rated speed 




(c)%10 rated torque, %10 rated speed 

Fig.8 20 hp induction motor estimated speed using reactive power MRAS scheme, 
applied torque and stator q-axis current 

V. EXPERIMENTAL RESULTS AND CONCLUSION 

A state estimator made by using MRAS has been tested experimentally as well. The 
experimental data; the real time stator voltages and currents are obtained from the setup is processed 
by Matlab in the host computer where the associated MRAS program is running. The outputs of the 
processing are displayed in Figs. 9. Gains of PI can be changed to improve the settling time, 
overshoot, rise-time, etc of the speed waveform while the system is going through the transient-state. 
The steady-state accuracy of MRAS meet the expectations and quite successful. Also, Fig.6.18 
shows the speed tracking performance of the back emf MRAS scheme. It is seen that this tracking 
performance of the speed estimator seems to be quite satisfactory. The simulations and experimental 
works show the great promise of the studied MRAS schemes 
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